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sive lepton-pairs in hadron reactions at fixed target and collider energies 
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within the context of next-to-leading order perturbative quantum chromo- 
dynamics. For values of the transverse momentum Qt greater than roughly 



><J ■ half the pair mass Q, Qt > Q/2, we show that the differential cross section 



is dominated by subprocesses initiated by incident gluons. Massive lepton- 
pair differential cross sections are an advantageous source of constraints on 
the gluon density, free from the experimental and theoretical complications of 
photon isolation that beset studies of prompt photon production. We com- 
pare calculations with data and provide predictions for the differential cross 
section as a function of Qt in proton-antiproton reactions at center-of-mass 
energies of 1.8 TeV, and in proton- nucleon reactions at fixed target and LHC 
energies. 
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I. INTRODUCTION AND MOTIVATION 



Both massive lepton-pair production, hi + h 2 — > 7* + X; 7* — > //, and prompt real pho- 
ton production, h\ + h 2 — > 7 + X are valuable probes of short- distance behavior in hadron 
reactions. In addition to the opportunities they offer for tests of perturbative quantum 
chromodynamics (QCD), the two reactions supply critical information on parton momentum 
densities. Massive lepton-pair production, commonly referred to as the Drell-Yan process M , 
provided early confirmation of three colors and of the size of next-to-leading contributions 
to the cross section differential in the pair mass Q. The mass and longitudinal momentum 
(or rapidity) dependences of the cross section (integrated over the transverse momentum Qt 
of the pair) serve as laboratory for measurement of the antiquark momentum density, com- 
plementary to deep-inelastic lepton scattering from which one gains information of the sum 
of the quark and antiquark densities. Inclusive prompt real photon production is a source 
of essential information on the gluon momentum density. At lowest order in perturbation 
theory, the reaction is dominated at large values of the transverse momentum px of the 
produced photon by the "Compton" subprocess, q + g — > 7 + q. This dominance is preserved 
at higher orders, indicating that the experimental inclusive cross section differential in pt 
may be used to determine the density of gluons in the initial hadrons [@-|5|]. 

There are notable similarities in the theoretical analyses of the two processes. At first- 
order in the strong coupling strength, a s , the Compton subprocess and the annihilation 
subprocess q + q — > 7 + g supply the transverse momentum of the directly produced prompt 
photons. Identical subprocesses, with the real 7 replaced by a virtual 7*, are responsible 
at 0(ct s ) for the transverse momentum of massive lepton-pairs. An important distinction, 
however, is that fragmentation subprocesses play a very important role in prompt real photon 
production at colliders energies. In these long-distance fragmentation subprocesses, the 
photon emerges from the fragmentation of a final parton, e.g., q + g — > q + g, followed by 
q -> 7 + X. 

Another significant distinction is that interest in hi + h 2 — > 7* + X has been drawn most 
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often to the domain in which the pair mass Q is relatively large, justifying a perturbative 
treatment based on a small value of a s (Q) and the neglect of inverse-power high- twist 
contributions (except near the edges of phase space). On the other hand, interest in prompt 
real photon production is directed to the region of large values of pr with small a s {pr)- 
Interest in the transverse momentum Qt dependence of the massive lepton-pair production 
cross section has tended to be focussed on very small values of Qt where the cross section is 
largest. Fixed-order perturbation theory |§ is applicable for large Qt, but it is inadequate 
at small Qt, and all-orders resummation methods H7HTTH have been developed to address 
the region Qt « Q- 

In this paper we focus on the Qt distribution for hi + h 2 — > 7* + X in the region 
where Qt is greater than roughly half of the mass of the pair, Q T > Q/2. We present 
and discuss calculations carried out in next-to-leading order QCD at both fixed target and 
collider energies. We show that the differential cross section in this region is dominated 
by subprocesses initiated by incident gluons. Massive lepton-pair differential cross sections 
are therefore a valuable, heretofore perhaps overlooked, source of constraints on the gluon 
density. We compare calculations with data and provide predictions for the differential 
cross section as a function of Qt in proton-antiproton reactions at a center-of-mass energy 
of 1.8 TeV, in proton-nucleon reactions at fixed-target momenta and proton-proton reactions 
at LHC energies. 

The requirement that both Q be large and Qt > Q/2 may well be punishing in terms 
of experimental rate, especially at fixed-target momenta. However, as long as Qt is large, 
the perturbative requirement of small « s (Qt) can be satisfied without a large value of Q. 
We therefore explore and advocate the potential advantages of studies of d 2 a / 'cIQcIQt as a 
function of Qt for modest values of Q, Q ~ 2GeV, below the range of the traditional Drell- 
Yan region. To be sure, one must avoid the region in Q of prominent resonance formation, 
e.g., the regions of the p, p', and J/ip resonances, among others. In addition, there are 
various backgrounds with which to contend at small Q such as the contributions to the 
event rate from prompt decays of heavy flavors, e.g., h\ + h 2 — > c + c + X; c — > I + X. These 
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heavy flavor contributions may be estimated by direct computation [12| and/or bounded 
through experimental measurement of the like-sign-lepton distributions. Previous studies of 
the connection between massive lepton-pair production and prompt photon production are 



reported in Refs. |0| and ]n|. 



As our calculations demonstrate, the Qt distribution of massive lepton pair production 
presents a valuable additional method for direct measurement of the gluon momentum dis- 
tribution. This method makes good use of a significant data base at collider and fixed target 
energies. The method is similar in principle to the approach based on prompt photon pro- 
duction, but it avoids the experimental and theoretical complications of photon isolation 
that beset studies of prompt photon production. 

In e + e~ and in hadron-hadron reactions at collider energies, prompt photons are observed 
and their cross sections are measured only if the photons are relatively isolated in phase 
space. Isolation is required to reduce various hadronic backgrounds including those from 
the electromagnetic decay of mesons, e.g., tt° — ► 27. The essence of isolation is that a cone 
of half-angle 6 is drawn about the direction of the photon's momentum, and the isolated 
cross section is defined for photons accompanied by less than a specified amount of hadronic 
energy in the cone, e.g., E™ ne < E milx = e^E^; E n denotes the energy of the photon. Instead 



of 5, the Fermilab collider groups use the variable R = y (A0) 2 + (A77) 2 , where Arj and A0 
denote differences of rapidity and azimuthal angle variables. Theoretical predictions will 
therefore depend upon the additional parameters and 8 (or R). The isolated cross section 
is not an inclusive cross section, and the usual factorization theorems for inclusive cross 
sections need not apply ]n|. Isolation removes backgrounds, but it also reduces the signal. 



For example, it reduces the contribution from processes in which the photon emerges from 
the long-distance fragmentation of quarks and gluons, themselves produced in short- distance 
hard collisions. 

The necessity to invoke phenomenological fragmentation functions and the existing in- 
frared ambiguity of the isolated cross section in next-to-leading order leaves some question 
about the extent to which isolated prompt photon data may be used for fully quantitative 



determinations of the gluon density. It is desirable to investigate other physical processes 
for extraction of the gluon density that are free from these systematic uncertainties. 

Fortunately, no isolation would seem necessary in the case of virtual photon production 
(and subsequent decay into a pair of muons) in typical collider or fixed target experiments. 
Muons are observed only after they have penetrated a substantial hadron absorber. Thus, 
any hadrons within a typical cone about the direction of the 7* will have been stopped, and 
the massive lepton-pair signal will be entirely inclusive. For this reason, we claim that the 
cross section for massive lepton-pair production at large Qt is an advantageous source of 
information on the gluon density. 

In Sec. II, we review the formalism of the next-to-leading order calculations of the trans- 
verse momentum distribution for both prompt photon production and massive lepton-pair 
production, demonstrating the analytic similarities and differences. In Sec. Ill, we present 
numerical calculations that establish the regions of Qt in which the gluon initiated Compton 
subprocess dominates the differential cross section d 2 a / 'alQalQx- Comparisons with data and 
predictions are provided in Sec. IV. Fixed-target data |TB[ are limited to the region Qt < Q, 
but they are useful for establishing the rough value in Qt — Q/2 above which resumma- 
tion may be ignored. Collider data |T7| extend to larger values of Qt- Our conclusions are 
summarized in Sec. V. 



II. MASSIVE LEPTON PAIR PRODUCTION AND PROMPT PHOTON 
PRODUCTION AT NEXT-TO-LEADING ORDER 

In inclusive hadron interactions at collider energies, hi + h 2 — > 7*+X with 7* — > 11, lepton 
pair production proceeds through partonic hard-scattering processes involving initial-state 
light quarks q and gluons g. In lowest-order QCD, at 0(oP s ), the only partonic subprocess 
is q + q — > 7*. At 0(a s ), both q + q — * 7* + g and q + g — > 7* + q participate, with the recoil 
of the final parton balancing the transverse momentum of the lepton-pair. These processes 
are shown in Figs. 1(a) and 2(a). Calculations of the cross section at order 0(a 2 s ) involve 
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virtual gluon loop corrections to these 0(a a ) subprocesses (Figs. 1(b) and 2(b)) as well as 
real gluon radiation contributions from a wide range of 2 — > 3 parton subprocesses (of which 
some examples are shown in Figs. 1(c) and 2(c)). 

The physical cross section is obtained through the factorization theorem, 

= £ J dx 1 dx 2 fi i (x 1 , ^)f h2 (x 2 , a^)-^-(«, Q, Qt, y; a*?)- (i) 

It depends on the hadronic center-of-mass energy S and on the mass Q, the transverse 

momentum Qt, and the rapidity y of the virtual photon; ///is the factorization scale of the 

scattering process. The usual Mandelstam invariants in the partonic system are defined by 

s — G°i + P2) 2 , t = {pi — P-y*) 2 , and u = (j>2 — £V ) 2 > where p± and p 2 are the momenta 

of the initial state partons and p 7 » is the momentum of the virtual photon. The indices 

ij G {qq, qg} denote the initial parton channels whose contributions are added incoherently 

to yield the total physical cross section. The partonic cross section <t?- (s,Q,Qt,V', H 2 ) is 

obtained commonly from fixed-order QCD calculations through 

d 2 <x 7 * d 2 a 1 *' {a) d 2 a 7 *' (b) d 2 a r,[c) 

atau dtdu dtdu dtdu 

The tree, virtual loop, and real emission contributions are labeled (a), (b), and (c) as are 

the corresponding diagrams in Figs. 1 and 2. The parameter /1 is the renormalization scale. 



It is set equal to the factorization scale [if = yQ 2 + Q\ throughout this paper. 

The cross section for h\ + h 2 — *■ 11 + X, differential in the invariant mass of the lepton 
pair Q 2 as well as its transverse momentum and rapidity, is obtained from Eq. (]l|) by the 
relation 

d 3 rr lT ( rv \ d 2 a r 



dQ 2 dQ\dy \3ttQ 2 J dQ\dy 
where Q 2 = (pi +P1) 2 , and pi,Pf are the four-momenta of the two final leptons. The so- 
called Drell-Yan factor a em / (3nQ 2 ) is included in all the numerical results for virtual photon 
production presented in Sees. Ill and IV. 

As stated in Sec. I, the partonic subprocesses that contribute to lepton pair production 
and to prompt photon production at finite transverse momentum are identical except for the 



fact that the photon is off-shell in the case of massive lepton pair production. As expected, 
the 0(a s ) cross section formulae for virtual photon production reduce to those for real photon 
production in the limit Q 2 — > |13[ . It was also noted earlier that fragmentation processes, 
which are not present in the Drell-Yan case, are important for real photon production. 
Taking into account direct and fragmentation pieces, we express the inclusive cross section 
for prompt photon production hi + h 2 — > 7 + X as 



dpldy 



Yl J dx 1 dx 2 fl 1 (x 1 , f/ f )f h2 (x 2 , 11}) (4) 

dz sdcr^' f ra 9 m 



dtdu 

The transverse momentum of the real photon is Pt- The fragmentation function D2(z,/j, 2 f ) 
is the probability for emission of a photon with momentum fraction z from parton k at 
the fragmentation scale \j? F . The function is a non-perturbative quantity that cannot be 
predicted by QCD but instead requires input from experiment at some reference scale. Its 
evolution to greater values of the scale is specified by perturbation theory. In this paper, we 
set fip = fif = fi. 

In prompt photon production, fragmentation contributes at leading order 0(a s ). Al- 
though the hard subprocess matrix elements af? f ragm in the fragmentation case are 0(a 2 ), 
the fragmentation function has a logarithmic dependence on the scale n 2 F and is effectively 
of 0(1/ a s ). When the subprocess matrix elements are convoluted with the fragmentation 
function, the product is of 0(a s ). 

At next-to-leading order 0(a 2 ), one must calculate virtual loop corrections to the leading 
order processes (Figs. 1(b) and 2(b)) and the three-body processes q + g — > 7* + q + g, 

q + q -> 7* + 9 + 9, Q + <? -> 7* + <? + <?, 9 + 9 -> 7* + 9 + 9, Q + q -> 7* + q + q, 
q + q' — > 7* + q' + q, q + q' — > 7* + q' + q' , and q + q — > 7* + q' + q' (Figs. 1(c) and 2(c)). The 
2 — > 3 parton matrix elements are integrated in 4 — 2e dimensions over the singular regions 
of phase space to expose the infrared and collinear singularities as 1/e-poles. These poles are 
then canceled against the poles present in the virtual corrections. In general, this procedure 



8 



applies to virtual photons as well as real photons. If the limit Q 2 — > is taken carefully, 
and no overall powers of Q 2 are factored out, the results for the virtual photon cross section 
will again reduce to those for the real photon cross section, as they do in leading order. 

There is, however, one particular case that has to be treated separately. In a next-to- 
leading order calculation of prompt photon production, a singularity is produced whenever 
the photon becomes collinear to a final state quark in the integration of the matrix ele- 
ments over phase-space. This collinear singularity must be subtracted and absorbed into 
the fragmentation function for a quark to produce a photon. The subtraction is done at 
an arbitrary fragmentation scale /xp, and hence the direct photon cross section acquires a 
logarithmic dependence on the scale fip that is partially canceled when fragmentation con- 
tributions are included. Fragmentation contributions are not present in the virtual photon 
case since the photon is off-shell, and its invariant mass Q regulates the singularity. This 
finite mass regularization leads to a logarithmic dependence of the Drell-Yan cross section 
on the virtuality of the photon or the invariant mass of the observed lepton-pair Q. As long 
as Q is large enough, these logarithms are small so that no resummation/f actor ization is 
necessary for a reliable prediction of the cross section. The Drell-Yan cross section will not 
have a fragmentation contribution except in the limit of very small Q. 

As an example of the discussion in the paragraph above, we consider the subprocess 
qq —>■ q'q'^i'J*) for prompt photon production and for the Drell-Yan case. The relevant 
diagrams are shown in Fig. 3(a) and 3(b) for initial- and final-state photon radiation. Since 
we are interested mainly in the structure of the final state in the region where a photon 
becomes collinear to a quark, this subprocess is typical and any conclusions drawn from it 
can be applied to the more complicated processes such as qg scattering. In this discussion 
we concentrate on the final-state collinear regions, but we use the full three-body matrix 
elements for both Drell-Yan and prompt photons H integrated fully over phase space, 
before factorization of any singularities. The exact forms for the phase space integrals can 
be found in Ref. H for prompt photons and in Refs. flT9"|,|r8fl for the Drell-Yan case. For the 
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direct photon cross section sd 2 aj- dtrect / dtdu, the phase space integrals yield a singular term 



of the form 



1 T ^ -(t 2 + u 2 )(2s 2 + 2st + t 2 + 2su + 2tu + u 2 ) 1 
-2iCv 7- — = -j{s,t,u) 



s(t + u) 4 

whereas the corresponding logarithmic term in the Drell-Yan case is 
~s + Q 2 -s 2 + \ 



(5) 



In 



K 



DY 



s + Q 2 -s 2 -\ 

1 (3Q 2 u 2 (u-t){t + u-2s 2 ) 3(s - Q 2 )u 2 (u 2 - t 



(6) 



A 5 



Q 2 + s - s 2 



+ 



1 / u(2ss 2 — 2s\ + 2s 2 t — su + 4s 2 u — tu — 3u 2 



A 3 



u 



1 1 (3s s 2 \ 
+ \\sU~J +U ) + Q2 + s-s 2 \2 

+ (t < — > u) 

' s + Q 2 - s 2 + A 



3s 5u 2u 2 ^ 

y + - 



In 



s + Q 2 - s 2 - A 



g(s,t,u,Q 2 ). 



The quantities K~ and K DY represent products of color, electric charge, and phase space 



factors. The function A = y (t + u) 2 — 4Q 2 s 2 ; s, t, and u are defined just below Eq. ([[]), and 
s 2 = s + t + u-Q 2 . 

The functions f(s,t,u) and g(s,t,u,Q 2 ) differ, of course, due to the Q 2 - dependence in 
the Drell-Yan case. In the limit Q 2 — > 0, however, g(s,t,u, Q 2 ) — > f(s,t,u), as expected. 
The logarithm in Eq. (6) becomes infinite in this limit and thus corresponds to the pole 1/e. 
For the prompt photon case the pole term is subtracted in some factorization scheme by the 
addition of a counter term. In the MS scheme the counter term is schematically 



< (5 » * 



(7) 



where a^ q is the hard 2-body subprocess matrix element, and P 1 / q {z) is the usual splitting 
function to produce a photon from a quark. This procedure replaces the pole term with new 
terms depending on ln(/i^). 



10 



To summarize this Section, we have shown that the cross section formulae for virtual 
and real photon production in hadronic collisions are intimately related in leading order 
0(a s ) as well as in next-to-leading order 0(a 2 ). If the logarithmic terms in Q 2 are first 
isolated, setting Q 2 = in the rest of the Drell-Yan cross section leads to a recovery of the 
prompt photon cross section except for the ln(//|*) terms and numerically small finite terms 
associated with the factorization of the collinear photon singularity into the fragmentation 
function. The situation is completely analogous to hard photoproduction where the photon 
participates in the scattering in the initial state instead of the final state. For real photons, 
one encounters an initial-state singularity that is factorized into a photon structure function. 
For virtual photons, this singularity is replaced by a logarithmic dependence on the photon 
virtuality Q |20[ . 



III. DOMINANCE OF THE COMPTON SUBPROCESSES FOR Q T > Q/2 

Having identified in Sec. II the subprocesses that contribute in next-to- leading order to 
massive lepton-pair production, hi + h 2 — > 7* + X, and to prompt real photon production, 
hi + hi — > 7 + X, and having established the correspondence between the analytic forms of 
the contributions to these two reactions, we turn in this Section to explicit evaluations of 
the differential cross sections as functions of Qt and pt, respectively, at both fixed-target 
momenta and collider energies. We work in the MS renormalization scheme and set the 
renormalization and factorization scales equal. We employ the CTEQ4M set of parton 



densities [^TJ and a two- loop expression for the strong coupling strength a s (p,), with five 



flavors and appropriate threshold behavior at fi = m&. For the massive lepton-pair case, a s 



is evaluated at a hard scale \i = yQ 2 + Q T , and for prompt real photon production, \i = px- 
Beginning with p + N — ► 7* + X at pi a b = 800 GeV, we present the invariant inclusive 
cross section Ed 3 a/dp 3 as a function of Qt in Fig. 4. Shown in this figure are the qq and 
qg perturbative contributions to the cross section at leading order and at next-to-leading 
order. For this calculation, we integrate the invariant inclusive cross section over the scaled 
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longitudinal momentum range 0.1 < xp < 0.3 and over the mass interval 5 < Q < 6 GeV, 
and we divide by the bin width in xp. To represent the target nucleon N, we use parton 
densities equal to 50 % of the proton densities and 50 % of the neutron densities, having in 
mind a comparison with deuterium data from Fermilab experiment E772, to be discussed in 
Sec. IV. 

The results in Fig. 4 show that the leading-order and the next-to-leading order contri- 
butions are about equal in size for Qt > 1 GeV in both the qq and qg channels, meaning 
that the cross section through next-to-leading order is about twice that of the leading-order 
value. For Qt < 2 GeV, the qq contribution exceeds that of qg channel. However, for 
values of Qt > 2 GeV, the qg contribution becomes increasingly important. As shown in 
Fig. 5(a), the qg contribution accounts for about 80 % of the rate once Q T — Q- The results 
in Fig. 5(a) also demonstrate that subprocesses other than those initiated by the qq and qg 
initial channels are of negligible import. 

In Fig. 5(b), we display the fractional contributions to the cross section as a function 
of Qt for a larger value of Q: 11 < Q < 12 GeV. In this case, the fraction of the rate 
attributable to qg initiated subprocesses again increases with Q T . It becomes 80 % or 
greater for Q T ~ Q/2. 

For the calculations reported in Figs. 4 and 5(a,b), we chose values of Q in the traditional 
range for studies of massive lepton-pair production, viz., above the interval of the J/ip and 
ip' states and either below or above the interval of the T's. 

For Fig. 5(c), we select the interval 2.0 < Q < 3.0 GeV. In this region, one would be 
inclined to doubt the reliability of leading-twist perturbative descriptions of the cross section 
da/dQ, integrated over all Qt- However for values of Qt that are large enough, a perturbative 
description of the Qt dependence of d 2 a / dQdQT ought to be justified. The results presented 
in Fig. 5(c) demonstrate that, as at higher masses, the qg incident subprocesses dominate 
the cross section for Q T larger than a few GeV. 

We return to pN reactions at fixed-target momenta in Sec. IV and discuss comparisons 
of our cross sections with data, but we address first the case of p + p — > 7* + X at collider 
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energies. 

In Fig. 6, we display the theoretical invariant inclusive cross section at the energy of the 
FNAL collider energy as a function of Qt for four regions of Q. In addition to the regions 
of Q mentioned above, we select a region of relatively large Q, 30.0 < Q < 35.0 GeV, to 
extend the region of interest of our analysis. For these results at collider energies, we average 
over the rapidity interval -1.0 < y < 1.0. The fractional contributions from the qg and qq 
subprocesses up through next-to-leading order are shown in Fig. 7. Evident in Fig. 7 is that 
the qg subprocess is the most important subprocess as long as Qt > Q/2. The dominance 
of the qg subprocess diminishes somewhat with Q, dropping from over 80 % for the lowest 
values of Q to about 70 % at its maximum for Q ~ 30 GeV. In addition, for very large Qt, 
the significant luminosity associated with the valence dominated q density in pp reactions 
begins to raise the fraction of the cross section attributed to the qq subprocesses. 

The calculations presented in Figs. 4-7 show convincingly that data on the transverse 
momentum dependence of the cross section for massive lepton-pair production at both fixed- 
target and collider energies should be a very valuable independent source of information on 
the gluon density. It is instructive to compare the results at collider energies with those 
expected for prompt real photon production. In Fig. 8, we present the predicted differential 
cross section for prompt photon production in the extreme cases: (a) all contributions, with 
full fragmentation taken into consideration (solid line), and (b) no fragmentation contribu- 
tions included (dashed line). The breakdown of these two cross sections into their qq and qg 
components is presented in Fig. 9. Comparing the magnitudes of the prompt photon and 
massive lepton pair production cross sections in Figs. 6 and 8, we note that the inclusive 
prompt photon cross section is about a factor of 400 greater than the massive lepton-pair 
cross section integrated over the mass interval 2.0 < Q < 3.0 GeV. This factor is attributable 
in large measure to the factor a em / (3ttQ 2 ) associated with the decay of the virtual photon 
to At current integrated luminosities, prompt photons have been observed with val- 



ues of pt extending to 100 GeV and beyond It appears that it should be possible to 



examine massive lepton-pair cross sections in the same data sample out to Qt of 30 GeV 
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or more. As may be appreciated from a comparison of Figs. 7 and 9, dominance of the qg 
contribution in the massive lepton-pair case is as strong if not stronger than it is in the 
prompt photon case. The statistical limitation to Qt of 30 GeV or so in the current data 
sample means that the reach in x g i uon , the fractional light-cone momentum carried by the 
incident gluon, is limited presently to 2Q T /\/~S ~ 0.033 or so, about a factor of three less 
than that potentially accessible with prompt photons in the same sample of data. It is 
valuable to be able to investigate the gluon density in the region Xgi uon ~ 0.033, and less, 
with a process that has reduced experimental and theoretical systematic uncertainties from 
those of the prompt photon case. 



IV. COMPARISONS WITH DATA AND PREDICTIONS 

In this Section we compare our calculations with data on massive lepton-pair production 
at large values of Qt, try to establish the domain in Qt over which the next-to-leading 
order perturbative calculations should be reliable, and offer predictions in addition to those 
presented in Sec. III. 

In Fig. 10, we show the invariant inclusive cross section Ed 3 a/dp 3 as a function of Qt 
for p + p -> 7* + X, with 7* -> fj+fT, at y/S = 630 GeV, with 2m M < Q < 2.5 GeV. 
These results are averaged over the rapidity interval -1.7 < y < 1.7. We present the next-to- 
leading order perturbative cross section along with its two major components, the qg and qq 
contributions. The theoretical expectation is in good agreement with the data published by 
the CERN UAl collaboration [17]]. Dominance of the qg component is evident over a large 
interval in Qt, whereas the qq contribution begins to be felt at very large Qt- 



Data obtained by the Fermilab E772 collaboration from an 800 GeV proton beam 

incident on a deuterium target are shown in Fig. 11 along with theoretical calculations. The 
solid lines show the purely perturbative next-to-leading order expectation, and the dashed 
and dotted curves pertain to results obtained from all-orders soft-gluon resummation and 
its perturbative expansion at low Q T - 
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The leading-order 0(a s ) and next-to-leading order 0{a 2 s ) perturbative calculations of 
the Qt distribution necessarily fail for Qt « Q-, where the analytic form of the fixed-order 
cross section behaves as 

da/dQ 2 T ~ ot s /Ql {l + a x a s \n 2 (Q/Q T ) + a 2 a 2 s ln 4 (Q/Q T ) + ...}, (8) 

In Eq. (H), <y s (fi) is evaluated at the scale fi = Q, and the are coefficients that do not 
depend on Qt- The expansion parameter of the series in Eq. (||) is a s ln 2 (Q /Qt), not a s , and 
it can become large at small Qt- The logarithmic series, associated with soft-gluon emission, 
is a vestige of the singularities that cancel at each order in perturbation theory between the 
real gluon emission and virtual gluon exchange diagrams. An elaborate procedure exists 
to sum the logarithmic series, known as resummation [0-0. We employ the formulation 



published by Arnold and Kauffman [10], along with the phenomenological non-perturbative 



functions of Ladinsky and Yuan ||11|| . We compare the resummed results with the purely 
perturbative results in order to try to establish the value above which we can be reasonably 
confident that fixed-order perturbation theory is adequate on its own. 

In Fig. 11, the dashed curves show the theoretical cross section obtained from resumma- 
tion. They are applicable only in the region Qt « Q. This restricted region of applicability 
arises because the only part of the fixed-order perturbative cross sections that is incorporated 
in the all-orders resummation formalism is the part that diverges as Qt 2 when Q T -> 0. The 
dotted curves in Fig. 11 show what is termed the "asymptotic" fixed-order contribution. It 
is obtained from an expansion of the resummed cross section in a power series in a s (Q). 
Mathematically, it represents the divergent part of the fixed-order perturbative cross section 
asymptotically in the limit Qt — > 0, and, as the coincidence of the dotted and solid curves 
shows in Fig. 11, this representation is exact at very small Qt- In our case, the dotted 
curves are the asymptotic answer through 0(a 2 ). 

Whereas resummation is clearly needed at small Qt, a fixed-order perturbative calcu- 
lation, without resummation, should be adequate at large enough Qt- In addition to the 
singular parts incorporated in the resummation, fixed-order perturbative cross sections in- 
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elude remainders that are non-singular at small Qt but are important at large Qt- An issue 
addressed in the literature is how best to connect the region in which all-orders resummation 
is pertinent to the region in which the simpler fixed-order result is fine. In efforts to restore 
the full content of fixed-order cross sections in the limit of large Qt, proposals have been 
made whereby the non-singular parts are reintroduced through "matching" prescriptions. 
The results of one approach to matching are sketched in Fig. 11 as the dot-dashed lines. The 
resummed result at very small Qt is joined smoothly to the purely perturbative fixed-order 
result at large Qt- 

Our purpose is to establish where the data on massive lepton-pair production may be 
used to test fixed-order perturbative QCD and to provide new constraints on the gluon 
density. The region of small Qt and the matching region of intermediate Qt are complicated 
by some level of phenomenological ambiguity. Within the CSS resummation approach, 
phenomenological non-perturbative functions play a key role in fixing the shape of the Qt 
spectrum at very small Qt, and matching methods in the intermediate region are hardly 
unique. For the goals we have in mind, it would appear best to restrict attention to the 
region in Qt above the value at which the resummed result (dashed curves in Fig. 11) falls 
below the fixed-order perturbative expectation (solid curves). A rough rule-of-thumb based 
on the calculations shown in Fig. 11 is Qt > Q/2. The published data at fixed- 



target momenta shown in Figs. 11(a) and 11(b) do not extend into this domain. Their reach 
is exhausted at just about the value in Q T at which the fixed-order calculation should be 
reliable. 

To establish the utility of massive lepton-pair production at fixed-target momenta for ex- 
traction of the gluon density, it would be valuable to obtain data in the traditional Drell-Yan 
region above Q ~ 4.0 GeV but whose reach in Qt extends into the fixed-order perturbative 
domain, as defined above. The massive lepton-pair production cross section increases as Q 
is reduced. In Fig. 11(c) we show the predicted cross section as a function of Qt for 2.0 
< Q < 3.0 GeV. The region Qt > 3.0 GeV is inviting. It is above the region in which 
resummation need be considered, and as shown in Fig. 5(c), it is a region in which the qg 
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subprocess dominates the final cross section. 

In Fig. 12, we provide predictions of the the invariant inclusive cross section Ed 3 a/dp 3 
as a function of Qt for p + p — > 7* + X at the LHC energy = 14 TeV. The curves show 
the next-to-leading order perturbative predictions integrated over four different intervals of 
Q and averaged over the rapidity interval -1.0 < y < 1.0. In all four cases, the contribution 
from the qg channel dominates at the level of 80 % or greater over the range of Qt shown. 
The LHC cross sections are about an order of magnitude greater than those at = 1.8 TeV 
over the range of Qt shown. 

V. DISCUSSION AND CONCLUSIONS 

Prompt real photon production is a source of essential information on the gluon mo- 
mentum density. At lowest order in perturbation theory, the reaction is dominated at large 
values of the transverse momentum p T of the produced photon by the "Compton" subpro- 
cess, q + g — > 7 + q. This dominance by gluon initiated subprocesses is preserved at higher 
orders, indicating that the experimental inclusive cross section differential in p T may be used 
to determine the density of gluons in the initial hadrons. There are notable similarities in 
the theoretical analyses of the massive lepton-pair production, the Drell-Yan process and 
prompt photon production. At first-order in the strong coupling strength, a s , the Compton 
subprocess and the annihilation subprocess q + q — > 7 + g supply the transverse momen- 
tum of the photons in both cases. The parallels are maintained in next-to-leading order, as 
summarized in Sec. II. 

In this paper we focus on the Qt distribution for hi + h 2 — > 7* + X. We present and 
discuss calculations carried out in next-to-leading order QCD at both fixed target and collider 
energies. We show that the differential cross section in the region Q T > Q/2 is dominated by 
subprocesses initiated by incident gluons. Dominance of the qg contribution in the massive 
lepton-pair case is as strong if not stronger than it is in the prompt photon case. Massive 
lepton-pair differential cross sections are therefore an additional useful source of constraints 
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on the gluon density. We compare calculations with data and provide predictions for the 
differential cross section as a function of Qt in proton-antiproton reactions at a center-of- 
mass energy of 1.8 TeV, in proton-nucleon reactions at laboratory momentum 800 GeV, and 
in proton-proton reactions at LHC energies. 

As long Q T is large, the perturbative requirement of small cz s {Qt) can be satisfied without 
a large value of Q. We therefore explore and advocate the potential advantages of studies 
of d 2 a I 'dQdQr as a function of Qt for modest values of Q, Q ~ 2GeV, below the range of 
the traditional Drell-Yan region. 

As our calculations demonstrate, the Qt distribution of massive lepton pair production 
offers a valuable additional method for direct measurement of the gluon momentum distri- 
bution. This method makes good use of a significant data base at collider and fixed target 
energies. The method is similar in principle to the approach based on prompt photon pro- 
duction, but it avoids the experimental and theoretical complications of photon isolation 
that beset studies of prompt photon production. No isolation would seem necessary in the 
case of massive virtual photon production (and subsequent decay into a pair of muons) 
in typical collider or fixed target experiments. Muons are observed only after they have 
penetrated a substantial hadron absorber. Thus, hadrons within a typical cone about the 
direction of the 7* will have been stopped, and the massive lepton-pair signal will be entirely 
inclusive. For this reason, we claim that the cross section for massive lepton-pair production 
at large Qt is an advantageous source of information on the gluon density. 

All-orders resummation is known to be important for the description of the Qt distri- 
bution of massive lepton-pair production at small and modest values of g-r- in this paper, 
we compare the resummed results with the purely perturbative results in order to try to 
establish the value above which we can be reasonably confident that fixed-order perturba- 
tion theory is adequate on its own. For the goals we have in mind, it would appear best to 
restrict attention to the region in Qt above the value at which the resummed result falls 
below the fixed-order perturbative expectation. A rough rule-of-thumb based on our calcu- 
lations is Qt > Q/2. The published data [|l6l at fixed-target momenta do not extend into 
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this domain. Their reach is exhausted at just about the value in Qt at which the fixed-order 
calculation should be reliable. To establish the utility of massive lepton-pair production at 
fixed-target momenta for extraction of the gluon density, it would be valuable to obtain data 
in the traditional Drell-Yan region above Q ~ 4.0 GeV but whose reach in Q T extends into 
the fixed-order perturbative domain. The experimental rates at the pp collider energy of 
\fS =1.8 TeV should be adequate over a wide range in Qt- 
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FIG. 1. (a)Lowest-order Feynman diagrams for the direct process q+q — > 7+5. (b) Examples of 
virtual gluon loop diagrams, (c) Examples of next-to-leading order three-body final-state diagrams. 
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FIG. 2. As in Fig. 1, but for the subprocesses initiated by the q + g initial state. 
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FIG. 4. Lowest order 0(a s ) (dashed lines) and next-to-leading order 0{ct 2 s ) (solid lines) per- 
turbative calculations of the invariant inclusive cross section Ed 3 a /dp 3 as a function of Qt for 
pN — ► 7*X at pi a b = 800 GeV, in the MS scheme. Contributions from the qg and qq channels are 
shown separately. The results are integrated over the scaled longitudinal momentum interval 0.1 
< xf < 0.3 and over the interval 5.0 < Q < 6.0 GeV, and they are divided by the bin width in xf- 
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FIG. 5. Contributions from the various partonic subprocesses to the invariant inclusive cross 
section Ed 3 a/dp 3 as a function of Qt for pN — > 7*X at pi a b = 800 GeV. The cross section is 
integrated over the scaled longitudinal momentum interval 0.1 < xp < 0.3 and over the intervals 
(a) 5.0 < Q < 6.0 GeV, (b) 11.0 < Q < 12.0 GeV, and (c) 2.0 < Q < 3.0 GeV, and divided by 
the bin width in xp- The contributions are labeled by qg (solid), qq (dashed), gg (dotted), qq~2 
non-factorizable parts (dot-dashed), and qq (wide dots). The qq and qq~2 results nearly coincide. 



28 



ft 

T3 



10' 



10 



10 



10 



pp->y*XatVS = 1.8TeV 
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FIG. 6. Differential cross sections Ed 3 a /dp 3 as a function of Qt for pp — ► 7*X at \/5 = 1.8 
TeV integrated over the rapidity interval -1.0 < y < 1.0 and over the intervals 2.0 < Q < 3.0 GeV 
(solid), 5.0 < Q < 6.0 GeV (dashed), 11.0 < Q < 12.0 GeV (dotted), and 30.0 < Q < 35.0 GeV 
(dot-dashed). The results of the integration have been divided by the bin width in y. 
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FIG. 7. Contributions from the partonic subprocesses qg and qq to the invariant inclusive cross 
section Ed 3 a/dp 3 as a function of Qt for pp — > j*X at \fS = 1.8 TeV. The cross section is 
integrated over the rapidity interval -1.0 < y < 1.0 and over the intervals 2.0 < Q < 3.0 GeV 
(solid), 5.0 < Q < 6.0 GeV (dashed), 11.0 < Q < 12.0 GeV (dotted), and 30.0 < Q < 35.0 GeV 
(dot-dashed). The results of the integration have been divided by the bin width in y. 
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FIG. 8. Differential cross sections Ed 3 a/dp 3 as a function of pr for real photon production 
pp — ► at = 1.8 TeV in the MS-scheme for two cases: no fragmentation terms included 
(dashed) and the inclusive case with full fragmentation included (solid). The rapidity y = 0. 
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FIG. 9. Separate contributions from the qq and the qg channels to the differential cross sections 
Ed 3 a /dp 3 as a function of pr for real photon production pp — > jX at y/S = 1.8 TeV in the 
MS-scheme for two cases: no fragmentation terms included (dashed lines) and the inclusive case 
with full fragmentation included (solid lines). The rapidity y = 0. 
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FIG. 10. Invariant inclusive cross section Ed 3 a/dp 3 as a function of Qt for pp — ► 7* + X at 
= 630 GeV, with 2m M < Q < 2.5 GeV, and averaged over the rapidity interval -1.7 < y <1.7. 
The next-to-leading order perturbative cross section (solid) is shown along with its two major 
components, the qg (dashed) and qq (dotted) contributions. The data are from the CERN UA1 
collaboration. 
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FIG. 11. Invariant inclusive cross section Ed 3 a/dp 3 as a function of Qt for pN — ► 7* + X at 
Plab = 800 GeV, averaged over the scaled longitudinal momentum interval 0.1 < < 0.3. The 
next-to-leading order perturbative cross section (solid) is shown along with the all-orders resummed 
expectation (dashed), the asymptotic result (dotted), and a matched expression (dot-dashed). The 
data are from the Fermilab E772 collaboration, (a) 5.0 < Q < 6.0 GeV, (b) 11.0 < Q < 12.0 GeV, 
and (c) 2.0 < Q < 3.0 GeV. 
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FIG. 12. Invariant inclusive cross section Ed 3 a/dp 3 as a function of Qt for pp — ► 7* + X 
at v 7 ^ = 14 TeV averaged over the rapidity interval -1.0 < y <1.0. The four curves show the 
next-to-leading order perturbative predictions integrated over 4 different intervals of Q: 2m M < Q < 
2.5 GeV, 5.0 < Q < 6.0 GeV, 11.0 < Q < 12.0 GeV, and 30.0 < Q < 35.0 GeV. 
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